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While most calculations on the properties of the ferromagnetic semiconductor GaAs:Mn have 
focussed on isolated Mn substituting the Ga site (Mnc* a ), we investigate here whether alternate 
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lattice sites are favored and what the magnetic consequences of this might be. Under As-rich 
(Ga-poor) conditions prevalent at growth, we find that the formation energies are lower for Mnc a 



, over interstitial Mn (Mm). As the Fermi energy is shifted towards the valence band maximum 

via external p-doping, the formation energy of Mn^ is reduced relative to Mnc a . Furthermore, 

in ' 

0^ ■ under epitaxial growth conditions, the solubility of both substitutional and interstitial Mn are 

O , 

strongly enhanced over what is possible under bulk growth conditions. The high concentration of 
Mn attained under epitaxial growth of p-type material opens the possibility of Mn atoms forming 
small clusters. We consider various types of clusters, including the Coulomb-stabilized clusters 



involving two M.iiQ a and one Mm. While isolated Mn^ are hole killers (donors), and therefore 



a 

q ! destroy ferromagnetism, complexes such as (MnG a -Mnj-Mnc a ) are found to be more stable than 



complexes involving Mnc a -Mnc a -MnG a . The former complexes exhibit partial or total quenching 

X" 

of holes, yet Mn^ in these complexes provide a channel for a ferromagnetic arrangement of the spins 
on the two Mnc a within the complex. This suggests that ferromagnetism in Mn doped GaAs arises 
both from holes due to isolated Mnc a as well as from strongly Coulomb stabilized MnG a -Mnj-Mnc7 a 
clusters. 

PACS numbers: PACS number: 75.50.Pp,75.55.-i,71.55.Eq 
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I. INTRODUCTION 

The discussion of the physics that underlies room-temperature ferromagnetism in 

transition-metal doped semiconductors has largely focussed on substitutional geometries, 
e.g. Mnc a site in GaAs. Indeed, there is a well-established tradition that 3d impurities in 
II- V semiconductors are largely substitutional while in Si they are mostly interstitials 
^J. Modern first-principles total-energy calculations afford testing of this classic paradigm. 
Recent experiments [5j find that Mn atoms occupy both substitutional as well as interstitial 
positions in GaAs. There have been suggestions from recent theoretical work J(| suggested 
that primarily surface energetics will funnel Mn atoms in to interstitial sites from surface 
adatom positions. While Mnc a behaves as a hole-producing acceptor, at the interstitial site, 
Mnj behaves as an electron-producing donor. Since ferromagnetism is mediated by free- 
carriers, Mnj could modify the magnetic properties from the case where only substitutional 
Mn sites were occupied. 

Using density-functional theory as implemented within plane-wave pseudopotential total 
energy method, we consider here bulk and epitaxial growth conditions, investigating isolated 
defects (Mnca and Mnj) and their complexes. We find that the Mn impurity in GaAs is 
stable in both substitutional and interstitial geometries depending on (a) the Fermi energy 
(which can be changed via external doping), (b) chemical potentials during growth and 
(c) bulk versus epitaxial growth conditions. The origin of these dependences is as follows: 
(a) As the formation energy of impurities that are neutral with respect to the lattice site 
they occupy (e.g. Mn^J does not depend on the Fermi energy (e^), the formation energy 
of positively charged impurities (e.g. Mn. + ) decreases as ep is shifted towards the VBM. 
Hence, the difference in the formation energies between Mn^ + and Mn^ a decreases with 
p-type doping, resulting in increased solubility of interstitial Mn. (b) Substitution of Ga by 
Mn involves the removal of a Ga atom and the introduction of a Mn atom at the site vacated 
by Ga. Thus, substitution is generally enhanced under Ga-poor, Mn-rich growth conditions. 
On the other hand the formation energy of Mn at an interstitial site does not depend on the 
Ga chemical potential. Thus, one may stabilize substitutional (interstitial) doping using Ga- 
poor (Ga-rich) growth conditions, (c) Solid solubility can be controlled thermodynamically 
using epitaxial instead of bulk growth conditions |7|. The absence of a substrate under 
bulk growth conditions allows the growing solid as well as its possible disproportionation 
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products to attain their free-standing lattice geometry. This is the case when the growth 
takes place from the melt as in Bridgman growth. Then, if phase separation occurs, the 
precipitate will take up its most stable crystal structure, i.e. MnAs in the NiAs structure. 
In contrast, under thin-film epitaxial growth conditions (as in MBE, MOCVD) competing 
phases such as phase separated MnAs are forced to be coherent with the GaAs substrate. 
As zincblende MnAs strained on GaAs is less stable than the NiAs phase of MnAs, phase 
separation is more costly under coherent epitaxial conditions, and one expects j?J less phase 
separation, hence enhanced solubility. We find the following: 

(i) Substitutional Mn has two stable charge states: neutral (Mn^ a ) and negatively charged 
(Mng a ) charge state which have 1 and holes, respectively. The calculated acceptor tran- 
sition E(0/-) between these states occurs at E v +0.13 eV, in good agreement with the 
experimental value of E„+0.11 eV [8j. Here E„ corresponds to the valence band maximum 
of the host material. 

(ii) The interstitial sites that Mn can occupy have either tetrahedral (coordinated to 
four As or four Ga atoms) or hexagonal symmetry. We find that Mn at the tetrahedral 
interstitial site coordinated by As is more stable than that coordinated by Ga, and exhibits 
a single charge state Mn 2+ for all values of the Fermi energy. The (0/+) and (+/2+) donor 
transitions are found to lie inside the conduction band, so, isolated Mn^ produce electrons 
that will compensate the holes created by Mn^. 

(iii) Under bulk growth conditions, the formation energy per Mn of substitutional 
Mn is AH(Mn^ a )=0.91+/iGa-/UMn eV=0.17 -HAs-^Mn eV, whereas interstitial Mn has AH 
(Mn 2 2+ )=0.55-^ M n+2e F eV per Mn. Here €p is the Fermi energy measured with respect to 
the valence band maximum of the host material, and fi^s and [iMn are the chemical poten- 
tials of As and Mn respectively. If we use maximally As-rich growth conditions (/Ms = eV) 
and //Mn=AH(MnAs), then we find AH(Mn^ a )=0.91 eV. As one dopes the sample p-type 
and ep approaches the VBM (e_F=0), the energy difference between the formation energies 
of Mn^. a and Mn 2+ reduces to 0.38 eV. It could decrease even further if 6f penetrates the 
valence band with doping, or if the growth conditions are made less As-rich. The interstitial 
concentration is then expected to further increase. This is confirmed by recent experiments 
|t| which use Ga-rich conditions for growth. 

(iv) Under epitaxial growth conditions, the formation energies of both substitutional 
and interstitial Mn decrease by 0.74 eV/Mn, so their concentrations increase concomitantly 
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leading to the possibilities of clusters. There is strong coulomb interactions between the 
oppositely charged constituents involving two substitutional (Mn^a) and one interstitial 
(Mnj); the cluster MnG a -Mnj-Mnc a is thus strongly stabilized and found to be more stable 
under p-type conditions than clusters involving three Mnc a . Epitaxial growth conditions 
increases the solubility of such Mnca-Mnj-Mn^a clusters, with formation energy of -0.15 +tp 
eV per cluster for the Q=+l charge state under As-rich conditions and /iA/ n =AH(MnAs). 

(v) The presence of interstitial Mn in the MnG a -Mnj-MnGa cluster provides a channel for 
the spins on the two substitutional Mn to align ferromagnetically even when there are no 
free carriers present in the cluster. We therefore conclude that ferromagnetism in GaAs:Mn 
can arise both from holes induced by isolated substitutional Mn atoms discussed previously 
(Dietl) as well as from charge compensated substitutional-interstitial clusters. 



II. METHOD OF CALCULATION 



The formation energy for a defect comprising of atoms a in the charge state q was 
computed using the density functional supercell method using the expression [loT ] 

AHj' q (e F ,fi) = E(a)-E(0) + Y^n a fi a a + q(E v + e F ), (1) 

where E(a) and E(0) are the total energies of a supercell with and without the defect a 
respectively. n a denotes the number of atoms of defect a transferred in or out of the reservoir 
(equal to 1 for an atom removed and -1 for an atom added), while \i a a denotes their chemical 
potentials. 

Total energies: The total energies of the charged supercells were computed by compensat- 
ing any additional charge on the impurity atom by a uniform jellium background and have 
been corrected for interactions between charges in neighboring cells using the Makov and 
Payne correction For isolated defects we used both the monopole as well as quadrupole 
corrections, while for composite defects we have added only the monopole correction to the 
total energy assuming all the charge to be localized at a single point. We use the static 
dielectric constant of GaAs (12.4) 12]. The quadrupole moment of the isolated defects was 
calculated as the difference between the moments of the supercell with the charged defect 
and that with the neutral defect. 
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Transition energies: The defect transition energy e(q, q') is the value of the Fermi energy 
6f at which AH a ' q (eF)=AH a ' q ' (ep). The zero of the Fermi energy is chosen as the valence 
band maximum E v of the pure host at the T point. 

Chemical potential limits: As the reservoir supplying the atoms could be elemental solids, 
or compounds formed from the elements, we express \x a a as the sum of the energy of the 
element in its most stable structure , and an additional energy fi a i.e \i a a = /x* + fjL a . The 
required ranges of fi a are determined by \iQ a < 0; [iun < 0; [Ia s < (no precipitation of 
solid elements) and by the formation energies of GaAs, and MnAs. The allowed values of 
chemical potential are such that GaAs is stable, i.e. fiGa + ^As = Affy(GaAs), the latter 
being the formation energy of zinc-blende GaAs. Further, as Mn should not precipitate as 
MnAs, we restrict [iMn + l^As < AH f (MnAs), the formation energy of MnAs in its most 
stable (NiAs) structure. For epitaxial growth conditions, the formation energy of zinc- 
blende MnAs lattice-matched to GaAs is considered. In this case we calculate the epitaxial 
formation energy, Ai//(MnAs) ep j, forcing the in-plane lattice constant of MnAs to equal 
that of GaAs, while the out-of-plane lattice constant, c, is allowed to vary. For coherent 
epitaxial growth the condition that MnAs should not form during incorporation of Mn in 
GaAs becomes fiMn + I 1 As < Aiff (MnAs) epi . 

The energies E(a), E(0), AfZj(GaAs), Af//(MnAs), AHf(MnAs) ep i and fi a are calculated 
within the density functional formalism, through the momentum-space pseudopotential total 
energy representation jl^j], using ultrasoft pseudopotentials [l^. The GGA-PW91 version of 
the exchange-correlation functional and no correction for the band gap underestimation 
was made. The calculations were performed over a Monkhorst-Pack 4x4x4 k-point grid for 
64 Q| and 216 atom supercells of GaAs using VASP Changing the k-point mesh from 
2x2x2 to 4x4x4 changed the formation energies by ~ 20 meV. Larger 256 atom supercells 
with 1x1x2 k-points were used for the calculations with clusters to ensure a larger separation 
between clusters. We used a plane wave cutoff of 227.2 eV for these calculations. Increasing 
the cutoff to 300 eV, changed the formation energies by ~ 10 meV. As the lattice constant 
of the supercell was kept fixed at the GGA optimised value for GaAs of a=5.738 A [18j, the 
internal coordinates were optimised. Our calculated (experimental) formation energies are 
AH f (GaAs) = -0.74 (-0.74), AH f (MnAs) = -0.74 (-0.61) eV and AH f (MnAs) epi ~ eV. 
For elemental Mn we assume the nonmagnetic fee structure j^, while for elemental Ga, we 
assume the base-centered orthorhombic structure. 
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The charge corrected |ll| Mn^a (0/-) transition as well as the difference in formation 
energies between Mn^ a and Mn^ + are given in Table I for supercell sizes of 64 and 216 
atoms. We see that changing the supercell size from 64 to 216 atoms lowers the acceptor 
energy by 30-50 meV and stabilizes Mn^, a over Mn^ + by 50-150 meV. The charge correction 
increases the acceptor energy by 60-90 meV and stabilizes Mn^ a over Mn 4 2+ by 250-350 meV. 



III. RESULTS 



A. Isolated substitutional Mn on the Ga site of GaAs 

Fig. 1 describes the formation energy AH(Mn^ a ) of neutral substitutional Mn in GaAs 
as a function of the chemical potentials has and fiMn- The shaded areas denote chemical 
potentials that produce unwanted products: (i) When /j,As becomes greater than zero (the 
cohesive energy of solid As) we have precipitation of elemental As as shown on the left 
hand side of Fig. 1. (ii) In the opposite limit, when /j,As takes more negative values than 
the formation energy AH(GaAs), we have maximally As-poor conditions and the host itself 
becomes unstable, as shown on the right hand side of Fig. 1. (iii) The diagonal lines in the 
main body of Fig. 1 denote different values of [iMn- When the chemical potential of Mn 
becomes greater than zero (the cohesive energy of solid Mn), metallic Mn will precipitate 
as shown in the bottom right corner of Fig. 1. Conversely, (iv) when /ij\/ n becomes equal or 
larger than AH(MnAs)-/x^ s , we will precipitate a secondary phase of MnAs. Clearly, since 
^Ga+I^As= AH(GaAs) and fJ-Mn+^As < AH(MnAs), one can keep the latter inequality even 
for moderately negative values of fiMn, provided that has is adjusted. The lines in Fig. 1 
show that the lowest AH(Mn^ a ) value is 0.91 eV (circle at bottom left corner). This can be 
attained at /Lt^ s =0 (maximally As-rich); /i^f n =AH(MnAs). Alternatively, the same solubility 
can be attained for less rich-As conditions, but richer Mn conditions, e.g. for / u j 4< j =-0.5 eV 
and /iAf n =-0.24 eV. 

Having described in Fig. 1 the stability of the neutral substitutional, we next describe 
in Fig. 2 the stability of the charged substitutionals. Here we chose the chemical potentials 
/iAs=0, //Mn=AH(MnAs) (denoted by the circle in Fig. 1) and vary the Fermi energy. We 
see that for p-type conditions, the lowest energy charge state is Mn^ a , whereas for higher 
Fermi energy the stablest charge state is Mn^. Table I gives the (0/-) acceptor transition 
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energy calculated with various supercell sizes with and without charge correction. The most 
converged (0/-) transition energy calculated for the 216 atom cell and corrected for charge 
interactions is E^+O.IS eV, in good agreement with the measured value of E„+0.11 eV |8[. 
Fig. 2 shows that under epitaxial conditions (right y-axis), the formation energy of Mng a is 
lowered by 0.74 eV. 

We next describe the electronic structure of Mn^a. In Figs. 3(a) and (b) we show the 
Mn d projected partial density of states (PDOS) for two charge states of substitutional Mn. 
The main features can be understood as arising from the hybridization between the anion 
dangling bonds generated by a Ga vacancy and the d levels on the Mn ion placed at the 
vacant site 3}. The Mn d ion levels are split by the tetrahedral crystal field into t 2 (d) and 
e(d). Exchange interactions further split these levels into spin-up (f) and spin-down (J.) 
levels. The t 2 (d) levels on the Mn atom hybridize with the levels with the same symmetry 
on the As dangling bonds, while the e(d) levels have no other states available for significant 
coupling [3] . Because the location of the Mn ion d levels is below the dangling bond levels, 
after hybridization, the deeper bonding t 2 states have dominantly Mn d character (refered to 
as CFR: "crystal field resonance"), while the higher antibonding t% states have dominantly 
As p character (refered to as DBH : "dangling bond hybrid"). These interactions lead to 
the energy level diagram depicted schematically on the left-hand-side of Fig. 4 showing a 
fully occupied, Mn-localized up-spin CFR of t 2 and e symmetries. At a higher energy we 
have the up and down-spin DBH states with t 2 symmetry. Because of the location of the Ga 
vacancy states t 2 (p) between the exchange split t 2 (d) states on the Mn, a negative exchange 
splitting is induced as a result of hybridization on the DBH states [21[ with t BBH below 
^dbh- As a resu h, the neutral substitutional defect Mn^ a has the electron configuration 
[t^e^]cFR (tft^)oBH, with a total magnetic moment [i=4 fis, and a hole in the t^ DBH orbital. 
This configuration corresponds to the multiplet 5 T 2 as observed in electron paramagnetic 
resonance (EPR) experiments (3]. The partial occupancy of the negative exchange-split 
DBH states stabilizes the ferromagnetic state over the antiferromagnetic state j^ . 

B. Isolated interstitial Mn 

Mn interstitial can occupy a site with tetrahedral symmetry (coordinated by four As or 
four Ga atoms) or a site with hexagonal symmetry. We have calculated the total energies of 
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Mn at these positions in a 64 atom cell of GaAs, and the results for the tetrahedral interstitial 
sites are given in Table II. The tetrahedral interstitial Mnj(As) coordinated by four As atoms 
is more stable than the one coordinated by four Ga atoms, with the difference being 0.16, 
0.31 and 0.31 eV for charge states q=l,2 and 3. In contrast, the hexagonal interstitial has 
0.62 eV higher total energy than the most stable Mn 2+ (As). Experimentally, the presence 
of interstitial Mn was detected by an analysis of the EPR spectrum as well as by 
Rutherford back scattering p. The distinction between the two types of interstitial sites 
(Mn-next to As vs Mn-next to Ga) is difficult to determine experimentally and involved an 
analysis of the experimentally measured contact interaction in terms of the covalency of the 
Mn-X bond. This analysis suggested that Mnj(Ga) was more stable, while our total energy 
calculations suggest that Mn;(As) is more stable. 

The formation energy of various charge states of interstitial Mn is shown in Fig. 2 for 
/iAs=0 and fiMn— AH(MnAs) . We see that the stable charge state is Mn. + for the full 
range of Fermi level, with maximum stability at e^=0. To compare the relative stability of 
Mn 2+ at 6f=0 with substitutional MiiQ a , we show in the upper scale of Fig. 1 the difference 
AH(Mn 2+ ) - AH(Mn^) between the formation energies of interstitial and substitutional Mn. 
We see that substitutional Ga is stabler on the left hand side of the figure, i.e sufficiently 
As-rich, whereas interstitial Mn is stabler at the right hand side of the figure, i.e. sufficiently 
As-poor. The energy difference is 

A#(Mn 2+ ) - AH(Mn° Ga ) = 0.38 + fi As + 2e F 

For /^4 S =0, the substitutional Mn are stabler by 0.38 eV, while for moderately As-rich 
conditions, say /j,a s = -0.4 eV, both defects have comparable formation energies. 

These results are in agreement with recent experiments using liquid phase epitaxy to 
introduce Mn in GaAs. Experimentally a decrease in hole concentration is found as the Mn 
concentration is increased. Under the Ga-rich growth conditions used, As antisites are not 
expected to be the dominant source of the observed compensation. Hence the major source 
of compensation is believed to come from Mn^ as expected for Ga-rich conditions from Fig. 
1. 

We next examine the electronic structure of Mn at a tetrahedral interstitial site. When 
Mn occupies a tetrahedral interstitial position, five of the seven electrons occupy the t^e^ 
CFR levels, with the remaining two going into the down-spin levels. This is evident from 
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the PDOS for the doubly-ionized Mn^ + shown in Fig. 3(c), where Mn^ + is found to have the 
configuration [t^e^]cFR with a magnetic moment of fi=5fiB- The central panel of Fig. 4(a) 
shows schematically the levels of Mrij. As the (0/+) and (+/2+) transitions are calculated 
to lie inside the GaAs conduction band (Fig. 2), we conclude that Mn^ produce free electrons 
in GaAs. 

C. Clusters of substitutional Mn 

Having dealt with the isolated limit, we investigated whether Mn atoms show a tendency 
to cluster. Recent experiments |24| on dilute magnetic semiconductors have found a strong 
tendency of the doped transition metal atoms to cluster and there has been some theoretical 
work 2j to support such observations. We consider As-centred clusters [(As)Mn n Ga4_ n ] 
with n=0, 1,2,3 and 4. 

Fig. 5 shows the formation energy of clusters made of three substitutional Mn atoms 
(S-S-S) at lattice locations (0,0,0), (a/2, a/2,0) and (0, a/2, a/2) in the 64 atom supercell. 
This corresponds to the n=3 cluster. Here a is the cubic lattice constant of GaAs. We 
see that the neutral cluster (3 Mn^a) having three holes is stable under p-type conditions, 
whereas the charged cluster (3Mnc a ) _ with 2 holes is more stable above 6f = 0.15 eV. 
The energies of the complex with 3(MnG a )° is 2.2 + 3(fj,Ga - I^Mn) eV while that of three 
noninteracting Mn Ga in their lowest energy charge state is 2.71 + 3(fiGa - l^Mn) eV. For 
epitaxial conditions AH ep j=0.02 +3(/i<3 a - nun) eV/cluster. Thus, as the formation energy 
is very low, the tendency for the Mn atoms to cluster is strongly enhanced under epitaxial 
growth conditions. 

In order to obtain a measure of the tendency to cluster, we calculate the clustering energy. 
The "clustering energy" S(n) is defined as the energy difference between n substitutional 
Mn atoms surrounding an As site [(As)Mn n Ga 4 _ n ; < n < 4] and n isolated well-separated 
constituents. Thus, 5E(n) = [E(n) — E(0)] — n[E(l) — E(0)}, where E(n) is the total energy 
of the supercell with As-centered clusters of n Mn atoms. We find that SE(n)= -228, -482 
and -794 meV per cluster of n=2,3 and 4 Mn atoms for a 64 atom supercell. The clustering 
energy changed to -519, -1069 meV for clusters involving 2 and 3 Mn in a 256 atom supercell. 
These results indicate a strong tendency for the neutral substitutional Mn atoms to form 
clusters. 
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D. Neutral complexes of Mna and Mn^ 

We considered the defect complex formed between Mn Ga (S) and interstitial Mn^ (I) 
denoted as (S-I-S)^, where Q is the total charge of the complex. The geometry for the 
complex had Mn Ga at (0,0,0) and (a/2, a/2,0) and Mnj at (a/2,0,0) in the 64 atom supercell 
of GaAs. We see in Fig. 5 that S-I-S exists in two charge states: When the Fermi energy is 
below E„+0.1 eV we have the stable structure is (S-I-S) 1+ , whereas when tp is above it, the 
stable structure is the neutral (S-I-S) . Thus, the donor transition for the cluster is at E„+0.1 
eV. Fig. 5 also shows that for Fermi levels below E^+0.22 eV, the S-I-S complex is stabler 
than the S-S-S complex. As for the interaction energy of the components of the complex: 
the formation energy of the non-interacting neutral (Q = 0) components of the complex is 
2E(Mn Ga )+E(Mn°)=4.31-|-2/i Ga - eV per 3 impurities, while the formation energy of 

the interacting neutral complex is 1.41+2/i Ga -3//Mn eV. This represents a ~ 2.9 eV per 3 
impurities stabilization over the non-interacting, neutral defects. The energy of the neutral 
complex measured with respect to the stablest lattice site occupied by isolated Mn under 
particular experimental conditions is found to be -574 meV for eV, /ZM n =AH(MnAs) 

and ep=0 eV. Hence this complex is strongly stabilized. 

The reasons for the stability of the (S-I-S) complex can be appreciated from Fig. 
4(a). Upon bringing together 2Mn Ga with Mn°, one electron drops from the higher en- 
ergy t? level of Mn, to the lower energy t® BH level of each substitutional site, resulting 
in [t^e^]cFR(tft^)DBH configuration at each Mn Ga site (Fig. 4(b)) which corresponds to 
Mn Ga . These conclusions are evident from our calculated DOS of the S-I-S complex, pro- 
jected on the I and S sites shown in Fig. 6. We find that for both Q=0 (Fig. 6(a)) and 
Q=2 (Fig. 6(b)) the I site has the configuration [t^e^]cFR or This substitutional-to- 

interstitial charge transfer lowers the energy of the complex by twice the separation between 
t? FR level of Mn, and t? BH level of Mn Ga . Furthermore, it creates a favorable Coulomb 
attraction between the components S~-I 2+ -S - of the complex. This energetically favorable 
substitutional-interstitial association reaction then eliminates the holes that were present in 



isolated substitutional Mn Ga and could explain the puzzling observation 25] of the existence 
of a far lower concentration of holes than Mn in GaAs. Alternate explanations such as the 
presence of As antisites |2j| as well as the presence of Mn atoms connected to six As atoms 
(as in the NiAs structure) have been offered. However, samples have been prepared where 
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the concentration of As antisites is too low to explain the observed compensation of holes. 
' further, six-fold coordinated Mn atoms have not been observed in Mn doped GaAs samples 
26|. 

E. Ferromagnetism of the (S — I — S)° complex 

The neutral complex has two Mn^ a and one intervening "d*" interstitial. We find that a 
ferromagnetic arrangement between Mnc a is favored in the complex (MnQ a —Mnf + —MnQ a )° 
by 176 meV. In contrast, our calculations for two Mng a atoms without the intervening 
interstitial atom finds that an antiferromagnetic arrangement of spins on the substitutional 
Mn atoms is favored by 108 meV. Thus Mnj is responsible for mediating a ferromagnetic 
interaction between Mn^. 

How does the presence of the interstitial Mn mediate the alignment of spins on the 
substitutional Mn? There are three possible arrangements for the spins on the Mn atoms 
making up the neutral complex - (S^S^)°, (S^S^) and (S^S^) . From our total energy 
calculations we find that the energy for the configurations (S^I^S^) and (S^S^) are higher 
by 563 meV and 176 meV, respectively, than the energy, E , of the ground state (SVS^) . 
(The energies changed marginally to 602 meV and 192 meV respectively when we increased 
the supercell size to 256 atoms.) The stabilization of the S^I^S^ magnetic arrangement can 
be understood using simple arguments: In the configuration (S^S^)°, as one spin channel 
is completely filled, there is no channel of hopping available for the electrons to delocalize 
and lower their energy. Thus, this is a high energy spin configuration with energy E +563 
meV. In contrast, in the configuration (S^S^)°, two channels of hopping are present; the 
first between the electrons on and S^, and the second between those on and P. This 
configuration is found to have the energy, Eo+176 meV. Likewise, the configuration (S^S^) 
which has energy Eo has two channels of hopping present between S and I. The dominant 
factor in determining the configuration which has the lowest energy are the hopping matrix 
elements - Vs,i between S and I and Vs,s between the two S. To a first approximation, these 
hopping matrix elements are determined by the separation between the atoms involved. As 
the distance between the two substitutional Mn atoms is v2 times the distance between S 
and I, the effective hopping matrix element between S is smaller. Hence the presence of an 
intervening Mnj provides a channel for the ferromagnetic arrangement of spins between two 
Mnc even in the neutral charge-compensated complex. In contrast, the presence of a closed 
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shell donor such as Asca between two Mn^a gives rise to an antiferromagnetic (or weakly 
ferromagnetic) interaction between Mn^a- 

How does the presence of the interstitial affect long-range ferromagnetism? In order to 
investigate this we introduced a hole-producing, isolated substitutional Mn atom at different 
lattice locations, indicated in Fig. 7, and investigated whether the spin on this isolated 
substitutional Mn atom prefers to align parallel or antiparallel with respect to the spins 
on the substitutional Mn atoms within the S-I-S cluster. We find that the substitutional 
Mn likes to align ferromagnetically with the Mn Ga of the cluster by 147, 214 and 81 meV, 
respectively for the positions 1, 2 and 3 (see Fig. 7). Hence the presence of the interstitial 
Mn forces a hole which is located ~ 12 A from the S-I-S cluster to align ferromagnetically 
and therefore contributes to the long-ranged ferromagnetism observed in these systems. 

As discussed earlier, the basic electronic structure of substitutional Mn in GaAs can be 
understood as arising from the hopping interaction between the Mn d states and the As p 
dangling bond states. Therefore, the coupling between two Mn atoms is through the As p 
states. It is strongest along the directions in which the "p — d n coupling of the Mn with the 
As states is the largest and decreases with distance along that direction. When there is a 
hole, the antibonding t 2 orbitals are partially occupied, and there is ferromagnetism. Hence 
in Fig. 7 spins on the Mn atoms at sites 2 and 3 prefer to align ferromagnetically with the 
spins on S when there is partial compensation. This could happen either because the S-I-S 
cluster is totally compensated, and a hole exists on the Mn at sites 2/3 or the S-I-S cluster 
is partially compensated. The mechanism stabilizing the ferromagnetic coupling between S 
and site 1 is the same as what we discussed for the (S-I-S) complex earlier and exists even 
when there is total compensation. 

F. Charged Mnca-Mrij-Mnca complexes: 

While the neutral complex has no holes, the Q—+1 and +2 complexes have 1 and 2 holes 
respectively (Fig. 4(c)) and a net magnetic moment of 4 and 3 /ib respectively. For Q=+2, 
Mn Ga adopts the configuration Mn^ a (Fig. 4(c)). We find that (Mn Ga -Mnj-Mn Ga ) 2+ prefers 
the ferromagnetic arrangement of spins on Mng„ by 286 meV, similar to the ferromagnetic 
preference (~ 305 meV) of Mn^ a -Mn^, a pair without an intervening Mn,. These results 
suggest the surprising fact that the spins on Mnc a align ferromagnetically in the charged 
complexes, almost as if Mn^ did not exist. The number of holes in the cluster is the same as 
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the number in the pair, though the number of Mn atoms are different. This is in agreement 

n 

with the experimental observation [27] where above a critical concentration of Mn, both the 
number of holes as well as the ferromagnetic transition temperature remain constant , while 
the magnetic moment per Mn atom decreases j^. The magnetic moments that we obtain 
for the Q=+l and +2 charge states translate into average moments of 1.33 \i& and 1 //# 
per Mn, while the uncompensated pair of Mnc a have a magnetic moment of 4 //g. In this 
regime where the T c is found to saturate, the average magnetic moment per Mn is found to 
vary from ~ 3 //# at a Mn concentration of 5.54 % to 1.74 at 8.3 %. 

Recent experiments [22( find that the T c of the as-grown samples increased after anneal- 
ing. This was interpreted as the migration of FM-reducing interstitial Mn to FM-enhancing 
substitutional positions. We investigated which clusters could break by annealing and pro- 
mote ferromagnetism. As the S-I-S complexes are rather strongly bound with respect to 
their constituents, we investigated instead complexes S-I, which are bound weakly (~ -196 
meV in the +1 charge state for fi^s, /UA/ n =AH(MnAs) and ej?=0 eV). We find antiferromag- 
netic spin arrangement in all Q=0, +1, +2 and +3 charge states considered. Thus, when 
these weakly bound S-I clusters are broken, depending on the charge state, there could be an 
increase in the number of holes and consequently the ferromagnetic transition temperature. 
On the other hand, S-I-S clusters appear to be stable and hence do not disintegrate under 
annealing. 
IV. SUMMARY 

Under As-rich conditions, Mn prefers to substitute the Ga site. As the growth conditions 
become less As-rich, or as extrinsic doping pushes ep towards and even below the VBM, 
the formation energy of interstitial Mn becomes competitive with that of substitutional 
Mn. Under coherent epitaxial growth conditions, when MnAs precipitates are forced to 
be coherent with the zinc-blende lattice, the formation energy of both substitutional and 
interstitial decrease. At this point, the solubility is large enough to form clusers. We find that 
S-I-S clusters are more stable than S-S-S clusters. S-I-S clusters are found to be strongly 
bound with respect to their constituents and exhibit partial or total hole compensation. 
While isolated Mnj behaves like a hole-killer and is expected to destroy ferromagnetism, in 
(Mnc a -Mnj-MnG a ) , the Mnj is found to mediate the ferromagnetic arrangement of spins on 
Miig a . The charged complex (Mnc a -Mnj-MnG a ) 2+ has a similar ferromagnetic stabilization 



13 



energy on the two Mn^a sites as in Mn^ a -Mng a cluster without Mn^ almost as if Mrii did 
not exist. Thus ferromagnetism in Mn doped GaAs arises from holes due to substitutional 
MiiGn, as well as from MnG a -Mnj-Mnc a complexes. 

This work was supported by the U.S. DOE, Office of Science, BES-DMS under contract 
no. DE-AC36-99-G010337. 
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FIG. 1: The formation energy of Miig a (left y-axis) as well as the difference in formation energies 
of Mn^ + and Mn^ a (top x-axis) are plotted as a function of fiAs (bottom x-axis) for different values 
of HMn- Here ep is fixed at the the VBM of the host. Regions where there is precipitation of the 
elemental solids as well as MnAs are also shown. 
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FIG. 2: The bulk (left y-axis) as well as epitaxial (right y-axis) formation energies for different 
charge states of isolated substitutional (S) and isolated interstitial (I) Mn calculated for a 64 
atom supercell under As-rich conditions. Acceptor transition for 216 atom supercell (Table I) is 
E^+O.IS eV. The chemical potentials are fixed at the points corresponding to the circle shown in 
Fig. 1. 
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FIG. 3: The £2 (upper panel) and e(lower panel) projected contributions to the Mn d projected 
partial DOS (a) for the q=0 and (b) -1 states of Mn^a as well as (c) the q=+2 charge state of Mnj. 



18 



(a) Non-interacting Mn° Ga and Mn ; 
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(b) Compensated, neutral Q=0 complex 
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(c) Charged, Q=+2 complex 
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FIG. 4: Schematic energy-level diagram for (a) neutral non-interacting substitutional (S) and 
interstitial (I) Mn impurities, (b) the compensated S-I-S complex and (c) the doubly charged S-I-S 
complex involving 2 substitutional and one interstitial Mn, where Q is the total charge of the 
complex. Open circles denote holes. 



Bulk and epitaxial formation of Mn clusters in GaAs 




e f (eV) 

FIG. 5: The bulk(left y-axis) as well as epitaxial (right y-axis) formation energies for different 
charge states of complexes involving two substitutional and one interstitial Mn compared with 
three substitutional Mn calculated for a 64 atom supercell under As-rich conditions. The chemical 
potentials are fixed at the points corresponding to the circle shown in Fig. 1. 
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DOS of Mn sites within the S-l-S cluster 




Energy (eV) 



FIG. 6: The ti (solid line) and e(dotted line) projected contributions to the Mn d projected partial 
density of states for Q=0 (left panels) and Q=+2 (right panels) of the complex projected onto Mn^ 
(top panels) and Mn<3 a (bottom panels). 
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FIG. 7: A single face of the 256 atom supercell of GaAs used in our calculations, where a is the 

cubic cell dimension. Positions 1,2 and 3 considered for the isolated Mnc a with respect to the 

cluster whose components are labelled S and I. 
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TABLE I: Acceptor transitions, Formation energies of Mnc a and Mnj for 64 and 216 atom super- 
cells of GaAs, with and without charge corrections 



Quantity 


64 atom cell with (without) 


216 atom cell with (without) 




charge correction (in eV) 


charge correction (in eV) 


Mn Ga (0/-) 


0.183 (0.094) 


0.133 (0.068) 


AH / (Mn^)-AH / (Mn^) 


0.382 (0.016) 


0.430 (0.17) 


/i M n=AH(MnAs),/XA s =0, e F =0 






AH / (Mn° Ga ) 


0.908 +flGa - fJ-Mn 


1.261 +jJLGa - I^Mn 



TABLE II: The formation energy for different charge states of isolated substitutional (Mnc a ) as 
well as interstitial Mn coordinated to four As atoms [Mnj (As)] or to four Ga atoms [Mnj(Ga)] 
where [i a denotes the chemical potential for atom a. 



Charge state 


Formation energy 




T d Mnj (As) 


T d Mnj(Ga) 


-1 


3.81-^Mn-CF 







2.45-/iMn 




+1 


l.l9-fi Mn +eF 


1.35-fiMn+^F 


+2 


0.18-flMn+^ F 


0A9-hmu +2e F 


+3 


0.24-/i M n+3e,F 


0.55-flMn+CF 
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